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Abstract 

We consider the oscillations of atmospheric neutrinos in the earth in the three- 
neutrino scheme with a Am^ hierarchy and a small admixture of the electron neutrino 
in the heavy mass eigenstate characterized by the mixing angle 6i^. We show that for 
Am^ ~ (0.5 — 3) X 10“^ eV^ indicated by the Super-Kamiokande data and sin^ 20 i 3 < 
0 . 2 , the oscillations of multi-GeV neutrinos in the subdominant i/g mode are 

enhanced by the MSW and parametric resonances. The parametric resonance, which 
occurs when the neutrinos cross the core of the earth, dominates for Am^ ~ (1 — 2 ) x 
10“^ eV^, sin^ 2013 ~ 0.06. The resonance matter effects lead to an observable excess 
of the e-like events with a specific zenith angle dependence even for small @ 13 . The 
up-down asymmetry of the multi-GeV e-like events can reach 15% for | cos0e| > 0.2 
and up to 30% for | cos0e| > 0.6, where 0e is the zenith angle of the electron. The 
resonance matter effects are relevant for the interpretation of the Super-Kamiokande 
data. 
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1 Introduction 


Neutrino oscillations in matter can be strongly enhanced if the matter density varies period¬ 
ically along the neutrino path. In this case the parametric resonance of neutrino oscillations 
can occur |]^, The probability of the transition of a neutrino from one flavor state to 
another may become close to unity even when the mixing angles (both in vacuum and in 
matter) are small. The parametric effects are further enhanced if the parametric resonance 
energy is close to that of the MSW resonance 0. 

The simplest realization of the periodic matter density distribution, which in addition is 
of practical importance, is the periodic step function (“castle wall”) prohle. In this case, the 
density modulation period L consists of two parts, and L^, which correspond to constant 
but different matter densities Nc and Nm- The evolution equation for neutrino oscillations 
in matter with such a density distribution allows for an exact analytic solution [Q. 

For the “castle wall” density prohle, the parametric resonance conditions are especially 
simple: the lowest-order (principal) resonance occurs when the oscillation phases <prn and 
(pc dehned through pi = 2TiLi/lm{Ni) [(i = m, c) with lm{Ni) being the oscillation length in 
the matter with density Np satisfy p|, ^ 

0c 0m ^ • (1) 

Recently it has been pointed out ^ that atmospheric neutrinos traversing the earth pass 
through layers of alternating density and can therefore undergo parametrically enhanced 
oscillations. Indeed, the earth consists of two main structures - the mantle and the core. 
The matter density changes rather slowly within the mantle and within the core but at 
their border it jumps sharply by about a factor of two. Therefore to a good approximation, 
one may consider the mantle and the core as structures of constant densities equal to the 
corresponding average densities (two-layer model) Q. Neutrinos coming to the detector from 
the lower hemisphere at zenith angles Qi, in the range dehned by cos0i, = (—1) (—0.837) 

comparison of neutrino oscillation probabilities calculated with such a simplified matter density 
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traverse the earth’s mantle, core and then again mantle. Therefore such neutrinos experience 
a periodic “castle wall” potential. Even though the neutrinos pass only through “one and 
a half” periods of density modulations (this would be exactly one and a half periods if the 
distances neutrinos travel in the mantle and core were equal), the parametric effects on 
neutrino oscillations in the earth can be quite strong i,i0- 

In 0 the effect of parametric resonance on possible ^sterile oscillations of atmo¬ 

spheric neutrinos was considered. The effect was found to be potentially important for zenith 
angle distributions of through-going and stopping muons produced by high energy neutrinos. 
The parametric enhancement of atmospheric neutrino oscillations <->■ was discovered 
numerically in [Q; however, in these papers the parametric nature of the enhancement was 
not recognized and possible consequences for atmospheric neutrino oscillations were not 
fully studied. 


The parametric enhancement can also take place for oscillations of solar neutrinos in 
the earth. It was realized in P, 0 that a sizeable enhancement of the regeneration effect 
for neutrinos crossing the core and having the energies lying between the MSW resonance 
energies in the core and in the mantle, found in 0, is due to the parametric resonance. 


In the present paper we discuss in detail possible manifestations of the resonance matter 
effects and in particular, of the parametric resonance in flavor oscillations of atmospheric 
neutrinos. The Super-Kamiokande (SK) collaboration has recently reported a strong evi¬ 
dence for neutrino oscillations in their atmospheric neutrino data [0, conhrming the pre¬ 
viously observed anomaly in the flavor composition of atmospheric neutrinos . 


The SK group has analyzed their data in the framework of two-flavor oscillations. They 
have shown that the data can be htted well assuming oscillations with the maximal 

or close to maximal mixing (sin^ 29 ~ 1) and Am^ = (0.5 — 6) x 10“^ eV^. Pure 

oscillations are practically excluded: they would have resulted in a signihcant zenith angle 
profile with those calculated with actual density profile provided by geophysical models shows a very good 
agreement (for recent discussions see, e.^., SD- 
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dependence of the e-like events in contradiction with the observations. In addition, the 
sizeable dehciency of muon-like events would require a strong (close to maximal) z/g — 


mixing, which is excluded by the CHOOZ experiment for most of the values of Am^ 
relevant for the atmospheric neutrino anomaly. At the same time, z/^ v^. oscillations 

with small mixing angles are still possible. Moreover, the data shows some excess of e- 
like events both in sub-GeV and multi-GeV samples and therefore is suggestive of the 
^ oscillations as a subdominant mode. In the data analysis by the SK collaboration 
the excess of e-like events was accounted for by up-scaling the overall normalization of 
the fluxes of atmospheric neutrinos by factors which are compatible with the uncertainty 
of the theoretical predictions. However, recent cosmic ray measurements by the BESS 
experiment indicate that the overall normalization of the atmospheric neutrino fluxes 
in the existing theoretical predictions was rather overestimated than underestimated. This 
makes the problem of the excess of e-like events more serious. In addition, the data, though 
not yet conclusive, seems to indicate some deviation of the zenith-angle dependence of the 
e-Iike events from the dependence that follows just from the zenith-angle dependence of the 
original flux in the absence of oscillations. Furthermore, it is difficult to explain the low 
value of the double ratio R = (e//i)/(e/p)Mc in a wide range of neutrino energies just in 
terms of the z/^ Vj- oscillations [Q. All this may indicate that the electron neutrinos are 
also involved in the oscillations. 


In this paper we study the atmospheric neutrino oscillations in the 3-flavor scheme with 
Vt being the dominant channel. We assume that the vacuum — v^. mixing angle 
satishes the upper bounds following from the GHOOZ data. We show that the MSW and 
parametric enhancements of atmospheric neutrino oscillations occur in the subdominant 
^ t'e niode, leading to observable effects despite the smallness of the mixing angle. 


There have been a number of studies of the atmospheric neutrino oscillations with the 
three neutrino mixing |^, [T^, [^, Q. However, their authors have either concentrated on 
gross characteristics of the atmospheric neutrino oscillations (such as the allowed values 
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of mixing angles and mass squared differences) [0, or considered only the cases when 
all mixing angles are large and/or is large [|^], or employed the constant-density 

approximation of the structure of the earth |^, or neglected the matter effects on neutrino 

In the hrst case the results are rather insensitive to the 


oscillations in the earth IS 


parametric effects, while in the last three cases these effects are missed altogether. 


Although the MSW enhancement effects are important and for a wide range of pa¬ 
rameters dominate the excess of the multi-GeV e-like events, we concentrate here on the 
parametric resonance effects for two reasons: (1) in contrast to the MSW enhancement 
which was widely discussed in the past the parametric effects in the oscillations of atmo¬ 
spheric neutrinos have not been studied in detail; (2) the parametric resonance modifies the 
MSW resonance peaks. 


The paper is organized as follows. In Sec. 2 we hnd the probabilities of oscillations in 
the three neutrino system in terms of the z/g — 1^3 oscillation probability which experiences 
the resonance enhancement. In Sec. 3 we consider the parametric resonance effects on the 
atmospheric neutrinos. In Sec. 4 we present the results of the numerical calculations for 
the zenith angle dependence and the up-down asymmetry of the e-like and /i—like events 
and confront the results with observations. In Sec. 5 we discuss our results as well as the 
prospects of observing the parametric resonance in the oscillations of atmospheric neutrinos. 
The details of the calculations of the cross-sections and the event rates are given in the 
Appendix. 


2 Three-flavor oscillations of atmospheric neutrinos 

We consider the three-flavor neutrino system with a hierarchy of mass squared differences 

Am2i -C Am32 ~ '^^31 • (2) 

We assume that Am ‘^2 = ~ 5 x 10“^ eV^ is relevant for the atmospheric neutrino 

oscillations whereas Am^i ~ 10“® eV^ allows one to solve the solar neutrino problem either 
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through the MSW effect or through the long range vacuum oscillations. 

The evolution of the neutrino vector of state Uf = (z/g, z/^)^ is described by the 

equation 


dt 






(3) 


where E is the neutrino energy and = diag{m\,rn\,rn\) is the diagonal matrix of 
neutrino mass squared eigenvalues. V = diagiVe, 0, 0) is the matrix of matter-induced 
neutrino potentials with Ve = V^GpNe, Gp and Nf. being the Fermi constant and the 
electron number density, respectively. The mixing matrix U, dehned through Uf = Uvrn 
where z/^ = (^i, ^ 2 , is the vector of neutrino mass eigenstates, can be parametrized as 


U = 1723t/i3f/l2. 


(4) 


The matrices Uij = Uij{ 6 ij) perform the rotation in the ij- plane by the angle 6 ij. We have 
neglected possible CP-violation effects in the lepton sector which are strongly suppressed in 
the case of the mass hierarchy 

Let us introduce new states z/ = (z/g, z> 2 , z/3)^ obtained by performing the U 23 - transfor¬ 
mation: Uf = f 723 Z 2 . The Hamiltonian H that describes the evolution of the z/ state can be 
obtained from (H) and (H): 

H = ^UM 2 M^uI,uI, + V . 


We get explicitly 


H 


s1^Aml2/‘^E + Ve 0 Si3Ci3Am32/2F^ 

0 cf2Am2i/2i? 0 

si3Ci3Am32/2£' 0 C33Am|2/2£' j 


(5) 


(ci 3 = COS 6 * 13 , Si3 = sin 6 * 13 , etc.) after the following approximations. Since Am‘^2 is in the 
range ^ 5 x 10“^ eV^ while m^i < 10“®eV^, the terms of the order S 32 Am 2 i/Am |2 were 
neglected. Also, Si 2 Ci 2 Am 2 i/ 2 ii^V'e 10 “^, so the ( 12 )-element in the matrix (^) (i.e. the 
mixing between the z/g and 1 ) 2 ) was neglected. 
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According to (H), the i >2 state decouples from the rest of the system and evolves inde¬ 
pendently. Therefore the S-matrix (the matrix of amplitudes) in the basis {ug, 1 ^ 2 , f's) has 
the following form : 




^ Age 0 AeS ^ 

0 A 22 0 

Ase 0 A 33 


( 6 ) 


where 


A 22 = exp(-i02), 


^12 


Am^iL 


2E 


(7) 


and L is the total distance traveled by the neutrinos. Notice that in our approximation 02 
does not depend on the matter density. The {ue, T’a) subsystem evolves according to the 2x2 
Hamiltonian [ue — z /3 submatrix in (^] determined by the potential I 4 , mixing angle ^13 and 
the mass squared difference Am‘^ 2 - L^t us denote by 


F 2 = |Ae3p = |A3gp = 1 - |Agg |2 = 1 - IH 33 P 


( 8 ) 


the probability of the Ue <->■ z /3 oscillations. As we will show in Sec. 3, it is in this channel 
that the oscillations are parametrically enhanced. 

The F-matrix in the flavor basis can be obtained from (|]) by U 23 rotation: U 23 SU 23 - 
Then the probabilities of flavor oscillations in the three neutrino system can be found as 
F(z/„ ^ U/ 3 ) = |(F235'Fj3)a/3p, which yields 


P{l^e 


^e) = 1 


2 , 


(9) 


F(z/e ^ Z/^) = F(z/^ ^ Z/g) = S23F2 , 

P{ye ^ ^r) = ^23-72 i 

P{j^t, ^ Z/^) = 1 - 5237^2 + 2 S 23 C 23 [Fe(e“*‘^M33) - 1 




\ 2 2 
^r) = ^23 C 23 


2 -P 2 - 2Re(e-*‘A23) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


The phase 02 is dehned in (^. The interpretation of the above results is straightforward. For 
instance, the z/^ <->■ z/g transition occurs via the projection of z/^ onto ^3 and 2 z/ - oscillations 
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z/g. The <->■ Ur transition can occur in two different ways: (i) projects onto U 2 , the 


latter propagates without transition and at the final time of evolution is projected onto Ur] 
(ii) Ufj, projects onto since z /3 oscillates into Ue, at the final time one should project the 
amplitude of the survival probability of onto Ur- The interference of these two transition 


amplitudes leads to the probability (|^. 


Using the probabilities d - (0) one can find the modifications of the atmospheric 
neutrino fluxes due to the oscillations. Let and be the electron and muon neutrino 
fluxes at the detector in the absence of oscillations. Then the fluxes in the presence of 
oscillations can be written as 


F — 


F = F^ 

^ e ^ e 


1-41- 


1 + P2{rs 


21 ' *23 


1 ) 


rslz, 


where 


r(^,0.) = 


P2 + 2S23C23 |^i?e(e *'^^^433) — 1 


(14) 

(15) 


is the ratio of the original muon and electron neutrino fluxes. 


It is interesting that one can have either an excess or a deficiency of e-like events de¬ 
pending on the values of r and S 23 . Indeed, the effect of oscillations on the electron neutrino 
flux is proportional to the factor (rs 23 ~ !)■ If one assumes r = 2 , there will be an excess of 
e-like events for 623 > 45° and a deficiency for 623 < 45°. The SK best fit was 623 = 45°; in 
this case there would be no deviation from the prediction for r = 2. However, for upward 
going neutrinos in the multi-GeV range r is typically 3 - 3.5 rather than 2 , so there should 
be an excess of e-like events even if 623 = 45°. In addition, notice that the SK analyses were 
performed in the two-flavor scheme, and the best-fit value of 623 may be somewhat different 
in the 3-flavor analysis. 

A final remark is that in the two neutrino mixing scenario, the probability of ^ Ur 
oscillations depends on sin^ 26 * 23 , and hence has an ambiguity 6*23 { 71/2 — 923 ). In fhe three 
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neutrino case the probability P(z/e <->■ depends on so that the study of an excess of 
the e-like events allows one to resolve the ambiguity. 


3 Parametric enhancement of neutrino oscillations 


The amplitude A33 and the probability P2 which enter into the expressions have 

to be found by solving the evolution equation for the (z/g, hs) system. The transitions in this 
system are the ones that undergo the resonance (parametric and MSW) enhancements. One 
can study properties of the resonance matter effects using the two-layer model of the earth’s 
density prohle. In the two-layer model P2 and A33 can be found explicitly in a compact 
form [0: 

/ : ( \ 2 


F2 = ( 2 sin sin 26 ^ Y + sin ^ sin 26 c j , 


( 16 ) 


Re[e *'^^33] = ^2 cos ^ Y — cos cos(<I) 


h) 


2 sin cos 26 m Y -|- sin ^ cos 26 c 1 sin($ — 02), 


( 17 ) 


where 

-VA _ 4’m 4^c . 0m . 0c /nn o/l \ o\ 

Y = cos cos — sm — sm — cos[ 26 m — 26 c) ■ (loj 

Here (pm and pc are the oscillation phases acquired by the neutrino system in the mantle 
(one layer) and in the core, respectively. They can be written as 


27rLi 

imm 


= Li AH (Vi), {i = m,c), 


( 19 ) 


where 


AH{Vi) = 


\ 


cos 2013 


Amo 


32 


2 E 


— Vi\ -I- sin 26 


Amo 


32 


13 ' 


2 E 


( 20 ) 


is the level splitting (difference between the eigenvalues of H), Vm and 10 being the potentials 
in the mantle and in the core. The angles 6m and 6c are the values of the z/g — ^3 mixing 







angle in matter of the mantle and the core respectively. They can be found from 


sin 26 i 


sin 26i3 


2EAH{Vi} ’ 


{i = m,c). 


( 21 ) 


The phase $ is given by the integral of (Am 32 / 4 i? + Ve/ 2 ) along the neutrino path between 
its production and detection points. In the two-layer model of the earth’s density prohle, it 
is given by 


$ 




Lm T 




2 


( 22 ) 


(we neglect the neutrino oscillations in the air which are unimportant for the range of 
the neutrino parameters of interest). For neutrinos crossing both the core and the mantle 
(sin^ 0 ,^ < = 0 . 299 ) the path lengths in the mantle and the core are determined 

by the relations 

Lm = R® cos Qy - \l{Rc/R®y - sin^ Qy ^ , (23) 

Lc = 2i?e ^(i?c/i?e)2 -sin^©,, . (24) 


Here i?® is the radius of the earth and Rc is the radius of the core. For sin^ 0 ,^ > 0.299 
neutrinos cross the mantle only and their path length is 2Lm = —2i?0 cos©,^. 


The physical picture of the oscillations and the resulting event rates depend crucially 
on the neutrino parameters Am‘^2 sin^ 26^13 as well as on the zenith angle 0 ,^. Let us 
consider the dependence of the oscillation probability on the neutrino energy. In the region 


E 


Ami 


= ( 0.8 - 4 ) • 10 ^' eV 




( 25 ) 


32 


neutrinos experience resonantly enhanced oscillations in matter. This interval is determined 
by the MSW resonance energies for oscillations in the core and in the mantle E^ and E^ 


Am|2 cos 26*13 


Er = 


2H 


(* = c, m) 


( 26 ) 


and by the resonance widths AE/Ef- ~ 2 tan 26*13. The MSW resonance enhancement 
leads to characteristic peaks in the energy dependence of the transition probability. The 
exact positions of the maxima of the peaks depend on the oscillation phases and in general 
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do not coincide with the resonance energies (^). Neutrinos having the trajectories with 
cos Qy > — 0.84 do not cross the core of the earth and therefore for such neutrinos only the 
MSW resonance enhancement of the oscillations in the mantle can occur. 

For cos Qiy < — 0.84 there is an interference between the oscillation effects in the core and 
in the mantle which strongly depends on sin 26^13. For sin^ 26^13 > 0.15 there is a signihcant 
overlap of the MSW resonances in the core and in the mantle. The interference leads to a 
rather complicated picture of neutrino oscillations in the overlap region with a modihcation 
of the MSW resonance peaks. The parametric resonance conditions are not fulhlled. 

For sin^20i3 < 0.15 the probability P2 as a function of E has three main peaks: two 
peaks with the maxima at ~ and ~ E^ due to the MSW resonance oscillations in the 
core and the mantle respectively, and a peak between them (fig. la). The latter is due to 
the parametric enhancement of oscillations. In what follows we shall call for brevity the 
resonance peaks due to the MSW effects for neutrinos oscillations in the earth’ mantle and 
core the mantle peak and the core peak respectively. 

The maximum of the parametric peak is at an energy Ep at which the resonance condi¬ 
tions (|I|) are satished; the analysis |P, of the parametric resonance condition shows that 
a significant parametric enhancement occurs only when 

^ (1 - 2 ) ■ IQi^ eV-', ( 27 ) 

i.e. Ep is indeed in the range E^ < Ep < E^. 

At the maximum of the parametric peak the transition probability (^) takes the value 

i 

Pr^ = sm^ 2 { 9 ,- 29 m) ( 28 ) 

provided that the resonance conditions are exactly fulhlled. Due to a small number of periods 
(“ 1.5 period”), the energy width of the parametric resonance AE/Ep is large. The transition 
probability decreases by a factor of two for [0] 

TT 

A0i = |0i - 7r| ~ - , (i = c,m). (29) 
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Thus, the resonance enhancement of neutrino oscillations can occur even for quite sizeable 
detuning of the phases 0 m,c- Numerically we get AE/Ep ~ 2 — 3 . 

Let us stress that the interference of the mantle and the core oscillation effects leads not 
only to the appearance of the parametric peak. It also modihes signihcantly the MSW peaks 
as compared with the peaks which would appear in the one layer cases without interference. 
In fact, the interference leads to a suppression of those peaks. 

The parametric resonance conditions (|^) constrain the allowed values of ^13. If these 
conditions are to be satisfied for all the neutrino trajectories that cross the core (including 
the vertical ones) one gets from (|lD the upper limit sin^ 26*13 < 7r^/(4(Lc)^„3. V^) ~ 0.04 0 . 
If the zenith angles close to 180 ° are excluded, the constraint becomes less stringent. For 
example, for sin^ > 0.12 one obtains sin^ 26*13 < 0 . 07 . 

The main features of the parametric peak are illustrated in fig. 1 . For cos 0 j, ~ — 1 
the parametric peak and the core peak partially overlap (actually the core peak appears 
as a shoulder on the low-energy slope of the parametric peak). With increasing cos 0 i, the 
parametric peak moves towards larger energies; it becomes well resolved from the MSW 
peaks in the interval cos 0 ,y = (— 0 . 94 ) (— 0 . 87 ) and eventually transforms into the mantle 

peak at cosQi, > — 0 . 85 . 

The relative strength (area) of the parametric and the MSW peaks depends on the value 
of 6*13. As follows from fig. lb, with decreasing ^13 the MSW peaks decrease faster than the 
parametric peak. This can be explained as follows. For a given value of ^13 and 0 j, there is in 
general some detuning at the parametric peak, i.e. the conditions (|l|) are only approximately 
satisfied. With decreasing sin^ 26*13 fhe maximal possible value (^) of the probability P2, 
which corresponds to the exact parametric resonance, decreases. At the same time, it turns 
out that the detuning becomes smaller which partially compensates the decrease in . 
Therefore in some range of ^13 the parametric peak lowers only moderately with decreasing 
sin^ 26*13. In contrast to this, the MSW peaks decrease rather quickly with sin^ 26*13. As a 
result, for sin^ 26*13 < 0.06 the parametric peak has the largest strength. 
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To summarize, the parametric peak is the most pronounced in the ranges of the param¬ 
eters characterized by eq. ( 0 ) and 

sin^ 2013 = (1 - 12) ■ 10"2 , cos 0^ = (-1) ^ (-0.84). (30) 


For the sub-GeV sample of e-like and /i-like events the relevant energies of neutrinos are 
E ^ 0.3 — 1.5 GeV. Then from (^) we hnd that for the sub-GeV events a signihcant effect 
of the parametric resonance is expected if the mass squared difference is in the range 


Am^2 ~ (1 - 10) ■ 10"^ eV^ . (31) 

The multi-GeV sample gets its main contribution from neutrinos with energies 1.3 —10 
GeV and the corresponding Am '^2 larger: 

Aml^ = (0.4 - 4) ■ 10-3 eV^ (32) 


Quite interestingly, the central value of this range coincides with the best £t value of Am| 


32 


which follows from the analysis of the contained events in SK [10 


The parametric resonance can play an important role for certain samples of events {e.g. 
multi-GeV) which pick up a relatively narrow neutrino energy interval. The number of events 
is determined by the integral over neutrino energy E of the oscillation probability folded 
in with the response function f{E): N oc J dEf{E)P{E). The response function describes 
the contribution of neutrinos with energy E to a given sample of events. In particular, the 
response function for the multi-GeV events in the Super-Kamiokande has a maximum at 
E ^ 2.5 GeV. The value of / increases rapidly with energy below the maximum but has a 
rather long tail above it; the energies at the half-height are E ~ 1.5 and 4 GeV. Thus the 
width of the response function is characterized by a factor of 2 - 3. 

As can be seen from £g. 1, for Am ^2 — (1-5 — 2) x IQ-^ eV^ and values of cos0i. 
varying in the range (—1) (—0.84) (i.e. covering the earth’s core), the parametric peak 

spans essentially the whole region of energies that corresponds to the peak of the response 
function of multi-GeV neutrinos, effect in the core). The mantle peak is in the tail of the 
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response function and therefore its contribution to the excess of e-like events is attenuated. 
For sin^ 2^13 < 0.06 the contribution of the parametric peak dominates for trajectories 
crossing the core of the earth. However, even for higher values of 013 (sin^ 2013 ~ 0.15) the 
contribution of the parametric resonance to the excess of e-like events can be comparable 
to that of the mantle and core peaks provided that Am 32 ~ (1 — 2) x 10“^ eV^. 

For Am 32 ^ 10“^ eV^ the overlap of the parametric peak and the peak of the response 
function is small and the main contribution to the excess of e-like events comes from the 
mantle peak. 

For Am |2 > (2 — 3) x 10“^ eV^, with increasing Am |2 the resonance effects weaken and 
the oscillation effects are essentially reduced to those of vacuum oscillations. 


4 Resonance effects and zenith angle dependence of 
e-like and /i-like events 

The resonance matter effects in the atmospheric neutrino oscillations can manifest them¬ 
selves as an enhanced excess of the e-like events with a specihc zenith angle dependence. 
We calculate the zenith angle distributions of the e-like and /i-like events for Am 32 = 
(0.1 — 10) X 10“^ eV^, sin^ 2023 ~ 0.7 indicated by the SK data and a wide range of values of 
sin^ 2013 taking the CHOOZ bounds into account. We examine whether the matter effects, 
and in particular the parametric resonance, can be relevant for understanding such features 
of the SK data as the excess of the e-like events and the asymmetries of e-like and /r-like 
events. 

The number of e-like or /x-like events with the detected charged lepton in the energy 
interval AEi and direction Qi (/ = e, fi) can be calculated as 

Jdn^i J dn, J dE,Ei{E,,n,)x 
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(Pa{E^,Ei,VLyi] 


- ((], © n,i)] e{Ei) , 


(33) 


dEi dElpi 

where Ei{Ey,Vty), (/ = e,/i) are the fluxes of neutrinos in the detector defined in eqs. (JT^, 
(H); El is the energy of the charged lepton. Up to the small geomagnetic effects the 
neutrino fluxes and therefore the charged lepton distributions depend only on the zenith 
angle: iV)(cos0;) = 2nNi{Qi). is the angle between the directions of the incoming 
neutrino and the outgoing lepton, (Pa/dEidQui is the neutrino charged current cross-section, 
e{Ei) is the charged lepton detection efficiency, and the integration goes over Ei > 1.3 GeV 
{El <1.3 GeV) for multi-GeV (sub-GeV) sample. 


In calculating the neutrino fluxes at the detector we have used the two-layer model of 
the earth’s structure [eqs. (|T^ ) - (Pl])]; the average densities of the core and the mantle 
were calculated for each neutrino trajectory using the actual density profile provided by the 
Stacey model [|^. For Amh < 10“^ eV^ the correction due to the phase 02 in eq. ([I5| ) is 
very small, and we have put 02 = 0 in the actual calculations. The details of the calculations 
of the cross-sections and the decay rates are described in the Appendix. 


The integration over the neutrino zenith angles and energies leads to a significant smear¬ 
ing of the Qi dependence. Indeed, the average angle between the neutrino and the outgoing 
charged lepton is about 15° —20° in the multi-GeV region and it is almost 60° in the sub-GeV 
range. Therefore the data do not give us direct information about the zenith angle depen¬ 
dence of Uf, and fluxes at the detector. In particular, a significant contribution to the 
vertical upward bin cos©; = (—1 -=—0.8) comes from the neutrinos which cross the mantle 
only and therefore do not experience the parametric enhancement of oscillations. Thus, the 
observed parametric enhancement effects are weakened as compared to the enhancement 
in the neutrino zenith angle distributions. This is especially true for the sub-GeV sample. 
Additional smearing of the neutrino zenith angle dependence due to the finite angular and 
energy resolution is relatively small. 
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In the antineutrino channels, matter suppresses oscillations | and consequently the para¬ 
metric effects are weak. Since neutrinos and antineutrinos of a given flavor are not distin¬ 
guished in the atmospheric neutrino experiments (he. only total v^. - 1 - v^. and fluxes 

are measured), the resonance enhancement effects are additionally diluted. 


Let us hrst consider the multi-GeV events. In Fig. 2 we show the zenith angle depen¬ 
dences of the e-like and /i-like events for a representative set of the oscillation parameters 
Am| 2 , sin^ 26^13 and sin^ 6 ' 23 . The e-like events exhibit an excess which hrst appears in the 
horizontal bin (cos0e — 0) and increases monotonically with decreasing cos0e. It reaches 
~ 20% in the vertical (upward) bin, cos0e — —1. In contrast to this, /x-like events exhibit 
a dehciency at cos 0 ^ < 0 . 


A convenient quantitative measure of the excess or dehciency of e-like and p-like events 
in the vertical bins is the up-down asymmetry 


<^^(& i ,& 2)=2 


N;^P{h,b2) - Nf°^^{b,,b2] 

K^{b,,b2) + Nf°^^{b,,b2] 


(I = e, fi) , 


(34) 


where 


iVr(&i,& 2 ) = f ''dcosQiNiiQi) , Nf°^^{b,,b 2 ) = dcos QiNi{Qi) , (35) 

J —^2 J bi 

and A)(cos 0;) are given in (^). The parametric as well as the MSW resonance enhancement 
ehects are largest in the vertical (upward) bins. In most of our calculations we use bi = 
0.6, 62 = 1 which corresponds to the SK binning. Notice that about 60% of neutrinos 
contributing to this bin have the trajectories which cross the mantle only and therefore do 
not undergo the parametric enhancement of oscillations. 


From the SK results 


101 we hnd 


A^/^(0.6,l) 0.22 ±0.21 . 


(36) 


The up-down asymmetry and the excess of the e-like events increase with increasing S 23 . In 
order to assess the maximal possible effects we therefore choose for most of our calculation 
the value S 23 = 0.75 which is close to the upper border of the values allowed by the SK data. 
^We assume Ato |2 > 0 throughout the paper. 
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In fig. 3a we show the dependence of the np-down asymmetry 74 ^/'^( 0 . 6 ,1) on Am ^2 
S 23 = 0.75 and different values of sin^ 26 * 13 . The dependence of on Aml 2 reflects the 
changing degree of overlap of the response function with the parametric peak as well as the 
MSW peaks (see Sec. 3). The maximum of the asymmetry (y4^/'^(0.6,1) ~ 0.15 — 0.25) 
is at Am |2 ~ (0.8 — 2) x 10“^ eV^. Thus, our calculations for sin^ 26*13 ~ 0.06 reproduce 
the central value of the up-down asymmetry observed by the SK (^). They also reproduce 
(within Icr) the experimentally observed excess of e-like events in the vertical bin. 

With decreasing sin^ 26*13 the asymmetry decreases and the maximum shifts to larger 
values of Am^. Indeed, for large 6*13 ( sin^ 26*13 rsj 0 . 1 ) the mantle peak gives the main 
contribution and the maximum of asymmetry at Am ‘^2 = 0-8 ^ 10 “^ eV^ reflects the position 
of this peak. With decreasing ^13 the parametric peak becomes relatively more important 
and the maximum of asymmetry shifts to Am 32 ~ (1.5 — 1.7) x 10“^ eV^ at sin^ 26*13 ~ 0.06, 
which corresponds to the position of the parametric peak. With further decrease of 6*13 the 
position of the maximum remains unchanged, in accordance with £g. lb. We hnd that for 
sin^ 26*13 0.025 the contribution of the parametric peak to the excess of the e-like events 

is about 60%. 

With decreasing Am \2 the asymmetry decreases since the mixing in matter (and con¬ 
sequently the oscillation probability) decreases. For Am ‘^2 < 0-3 x 10“^ eV^ it approaches 
the asymmetry due to the geomagnetic effect without oscillations. 

For Am ‘^2 > 3 x 10“^ eV^ the asymmetry decreases with increasing Am ‘^2 reasons, 

(i) The matter enhancement of mixing disappears and the oscillation effect is essentially 
reduced to that of vacuum oscillations governed by the small sin^ 26*13 (the oscillations due 
to the matter splitting between the levels of the two light eigenstates are suppressed by 
the factor sin^ 26*13 ■ sin^ 26 * 23/4 |^); (ii) for Am \2 > 5 x 10“^ eV^ the oscillations become 
important for down-going neutrinos too and the up-down asymmetry goes to zero. 

In fig. 3b we show the asymmetry of p-like events, which is opposite in sign to that for the 
e-like events. The absolute value of the asymmetry becomes maximal at Am ‘^2 ~ 0-3 x 10“^ 
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eV^. This corresponds to the situation when the average distance in the vertical bin (~ 10^ 
km) equals half of the vacuum oscillation length for a typical energy of multi-GeV neutrinos 
E ^ 2 GeV. The matter effects are important in the range Aml 2 ~ (0.5 — 5) x 10“^ eV^ 
where the absolute value of the asymmetry increases with sin^ 29i3 just as in £g. 3a. We 
show by dotted line the asymmetry which would be expected in the case of pure 
Ur oscillations with the same S23 = 0.75. (For 2 x 10“^ ^ Am| 2 /eV^ $ 10“^ it can be 
roughly estimated assuming that in this range the oscillations of upward going neutrinos 
are fully averaged whereas the downward going neutrinos do not oscillate at all, which yields 
^ 26*23/(2 — 0.5 sin^ 26 * 23 ) = —0.46.) The difference between the dotted line and the 

others shows the enhancement of the asymmetry due to an additional channel of oscillations 
^ iu the three neutrino system. In the case of •«-*• u^ oscillations with maximal 
mixing the asymmetry is larger: ~ —0.67. Both values are compatible with the SK 

result, —0.56 ± 0.15 |I0[ . 


As we have pointed out above, the parametric enhancement of neutrino oscillations 
occurs for neutrinos crossing the earth’s core (cosG^, < —0.84). Therefore the largest up- 
down asymmetry is achieved when the binning enhances the contribution of the core-crossing 
neutrinos. In £g. 3c we compare the electron asymmetries in the bins 0.84 < | cos0j,g| < 1, 
0.60 < I cos © 1,^1 < 0.84 and 0.60 < |cos0yj < 1. One can see that the asymmetry is 
largest for the hrst bin which has the maximal contribution from the neutrinos crossing the 
core of the earth. The maximum of asymmetry is achieved at Am |2 — (1.5 — 1.7) x 10“^ 
eV^ which is typical of the parametric resonance at multi-GeV energies. The position of the 
peak of the asymmetry in the second bin which is dominated by the mantle-only crossing 
neutrinos (Am 32 — 0.7 x 10“^ eV^) reflects the position of the MSW peak in the mantle. 
The parametric peak is higher than the MSW one because for small sin^ 2^13 the parametric 
effects dominate. In the third (largest) bin, 0.60 < | cos0j,^| < 0.84, the parametric effects 
are signihcantly diluted compared to the hrst bin. 


To summarize, the physical picture of the neutrino oscillations depends crucially on the 
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value of Am‘^ 2 - The whole range of Am ‘^2 can be divided into three parts: 

( 1 ) the region of the vacuum oscillations: Am \2 ~ 3 x 10“^ eV^ for multi-GeV neutrinos 
(here there is also a small effect due to the matter-induced level splitting of the light states); 

(2) the resonance region (two MSW resonances and the parametric resonance): Am \2 — 
(0.5 - 3) X 10-3 eV^; 

(3) the region of matter suppressed oscillations Am |2 ~ 0.5 x lO-^ eV^. 

For sub-GeV neutrinos the corresponding regions are shifted by about a factor of 3 to 
smaller values of Am\ 2 - 

In £g. 4a, we show the dependence of the up-down asymmetry y4^/^(0.6,1) on sin^ 2^13 
for s |3 = 0.75 and several values of Am‘^ 2 - The dashed curve corresponds to the value 
Am 32 = 1.7 X 10-3 eV^ from the resonance region. The asymmetry (and the excess) of the 
e-like events rapidly increases with sin^ 29i3 in the region sin^ 2^13 < 0.10 which corresponds 
to the parametric resonance, and then increases more slowly. 

The solid line corresponds to a value of Am |2 close to the vacuum oscillation region. 
For sin^ 2013 > 0.08 the asymmetry increases almost linearly with sin^ 20 i 3 , as is expected 
in the case of vacuum oscillations. Here, too, the GHOOZ bound is important: for the 
allowed values (sin^ 20 i 3 < 0.16) the asymmetry is smaller than that for the value Am ^2 — 
1.7 X 10-3 eV^ from the resonance region. The dot-dashed line represents the asymmetry in 
the region of the matter suppressed oscillations: the effective mixing is suppressed even for 
large sin^ 20 i 3 and the asymmetry is relatively small. In £g. 4b, we show the corresponding 
asymmetry in the /x-Iike events. 

As follows from ( 0 ). the asymmetry in the e-like events is proportional to the factor 
(fs^s — 1), where r ~ 2.5 is the ratio of the original muon and electron neutrino fluxes 
averaged over the zenith angles and energies in the multi-GeV sample. According to this, 
the asymmetry due to the oscillations increases almost linearly with S 23 (hg- 5a); it becomes 
zero at S 23 ~ r~^ ~ 0.4, where the total asymmetry, A ~ 0.05, equals the small asymmetry 
due to the geomagnetic effects. 
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In fig. 5b, we show the dependence of the up-down asymmetry for p-like events on 
S 23 . For pure <->• Vr vacuum oscillations (dotted curve), the up-down asymmetry is a 
symmetric function with respect to S 23 (1 — " 523 ). Matter effects break this symmetry and 

the breaking increases with in accordance with our previous analysis. 

With increasing S 23 the excess of the e-like events gets enhanced. At the same time, for 
6^23 > 45°, the probability of the Vr oscillations decreases with increasing The 

additional channel of oscillations <->■ z/g does not compensate for this decrease, and so the 
asymmetry and the total suppression of the p - like events become smaller (see £g. 5 a, b). 

The zenith angle dependence of the multi-GeV e-like events with and without z/g z/^ 
oscillations is shown in fig. 6 , along with the SK data for 535 days. One can see that taking 
into account the z/g <->■ z/^ oscillations improves the £t of the data but cannot fully explain 
the excess of the e-like events unless the overall normalization of the atmospheric z/g and z/^ 
fluxes is increased. 

In fig. 7 the iso-asymmetry curves for multi-GeV e-like events are plotted in the (sin^ 2^13, 
Am 32 ) plane, along with the constraints from the GHOOZ experiment (shaded area). The 
behavior of the curves can be understood from the preceding discussion. The large-Am |2 
region of the plot corresponds to the vacuum oscillations; vertical and near-vertical lines are 
due to the averaged oscillations with no or little dependence of probability on Am 32 . Matter 
effects increase the asymmetry in the region Am |2 — (0.5 — 3) x 10“^ eV^: the iso-asymmetry 
curves are “pulled” towards the region of small values of sin^ 20 i 3 . For Am ‘^2 ~ x 10“^ 
eV^ matter suppresses the z/g z/^ oscillations. 

The parameter space of the resonance region can be divided into three parts: (i) Am \2 ~ 
10 “^ eV^, where the mantle resonance dominates; (ii) Am \2 ~ 1-5 x 10 “^ eV^, sin^ 2^13 < 
0.06, where the parametric resonance gives an important contribution; (hi) Am |2 ~ 1-5 x 
10 “^ eV^, sin^ 26^13 — 0.06 — 0 . 12 , where parametric peak gives a smaller contribution which, 
however, is comparable to that due to the mantle resonance. The core resonance gives a 
smaller effect. 
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For sin^ 2^13 > 0.15 there is a complex interference of the MSW resonance effects in the 
mantle and core. 

The CHOOZ bound excludes the part of the parameter space that corresponds to large 
asymmetry, 1) > 0.28. There is a local maximum of the asymmetry 74^/^(0.6,1) ~ 

0.27 at sin^ 2^13 ~ 0.14, Am ‘^2 — 10“^ eV^. Values close to this can also be achieved at the 
same Am \2 — 10“^ eV^ but large mixing angles sin^26'i3 ~ 0.5 — 0.9. From £g. 7 it follows 
that the asymmetry y4^/'^(0.6,1) can be as large as 0.22 in the range of parameters (ii), i.e. 
rather close to the maximal possible value indicated above. 

The parametric enhancement of the ^ oscillations makes it possible to have a 
sizeable asymmetry even for very small values of the mixing angle 013: the asymmetry can 
be as large as 0.15 even for sin^ 2^13 = 0.02 provided that Am ^2 li^s in the range (1—2) x 10“^ 
eV^. Notice that in the absence of matter effects one would expect the asymmetry to be an 
order of magnitude smaller. 

The iso-asymmetry plot of fig. 7 has been obtained for the hxed value of S23 (^23 = 0.75); 
the magnitudes of the asymmetry for other values of can be easily found using the relation 
= 2x1 {2 x), where x = — 1)^2- 

As we have pointed out above, the largest asymmetry allowed by the CHOOZ con¬ 
straints, A^/-^(0.6,l) ~ 0.28, is achieved in a small region around Am \2 — 0.8 x 10 ^ eV^ 
and sin^ 26*13 — 1 (see £g. 7). Future reactor experiments, and in particular KAMLAND 
||20|| , will be able to probe this range of parameters. If no oscillations are found, the values 
sin^ 2013 >0.1 will be excluded. In this case the largest possible asymmetry would corre¬ 
spond to the small-013 region where the parametric resonance effects play an important role. 
It should also be emphasized that in case of the negative result of KAMLAND the studies 
of the excess of the e-like events in atmospheric neutrinos would be a unique way to probe 
the parameter range sin^ 20i3 < 0.1, Am ^2 > 5 x 10“^ eV^. 

The calculated event rates are sensitive to the response function used. The latter de¬ 
pends on the event selection criteria, detection efficiency and other features of detector. 
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etc.. The results given above correspond to the response function described in Sec. 3. 
The response function with the maximum (and median energy) shifted by about 25% to 
higher energies compared to what we used would shift the iso-asymmetry contours by about 
25% to larger Am^. Notice that in this case the region of large asymmetries at large mix¬ 
ing angles will be completely excluded and the largest asymmetry would be achieved at 
sin2 20i3 ~ 0.15. 

Let us now consider possible effects in the sub-GeV sample. In fig. 8 we show the zenith 
angle dependence of the sub-GeV events for Am32 = 0.3 x 10“^ eV^ which corresponds to 
the maximum of the resonance effects. Due to the strong smearing, the excess has a rather 
weak zenith angle dependence. Moreover, even in the vertical bin it does not exceed 10%. 
The up-down asymmetry is smaller than 0.03. Notice that in the sub-GeV region f is closer 
to 2 than in the multi-GeV region and the factor fs23 — 1 in (0) leads to an additional 
suppression of the transition probability. 

In fig. 9a we show the dependence of ratio N^/N^ of the total numbers of the sub- 
GeV events on Am32 with and without oscillations. The matter oscillations enhance the 
excess of e-like events in the range Am ‘^2 = (0.1 — 1) x 10“^ eV^ and the maximum is at 

Am‘^2 — 0-3 X 10“^ eV^. 

For Am ‘^2 > 10“^ eV^, which corresponds to the maximal effect in the multi-GeV sample, 
the oscillation effect in the sub-GeV sample approaches that of the vacuum oscillations. For 
sin^ 26*13 < 0.1, we find that the total excess of the e-like events is below 5% and the 
asymmetry is below 3%. 

The ratio N^/N^ is mainly determined by the —>■ oscillations. In the 2-flavor case 
(6*13 = 0) these oscillations are unaffected by matter. Therefore in the 3-flavor case for small 
values of 6*13 the matter effects on are relatively small (fig. 9b). 
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5 Discussion and conclusions 


The excess of e-like events both in multi-GeV and sub-GeV samples observed by the SK 
experiment may indicate that the electron neutrinos are involved in the oscillations of at¬ 
mospheric neutrinos. In this connection we have considered the oscillations of atmospheric 
neutrinos in the 3z/ scheme. Assuming the mass hierarchy Am ^2 ~ ^ have 

derived simple analytic expressions for the oscillation probabilities in the 3z/ system in terms 
of the oscillation amplitudes in the 2 u system (z/g, ^ 3 ), where z /3 is a linear combination of 
Ufj_ and Ur- For the amplitudes of the z/g z /3 oscillations in the earth analytic expressions 
obtained in the two-layer model of the earth’s structure were used. 

Let us summarize our main results. We have shown that the range of the neutrino 
parameters sin^ 26*13 ~ 0.2, Am \2 ~ (0.5 — 3) x 10“^ eV^ is the resonance range for multi- 
GeV events in which the oscillations in the subdominant z/g <->■ z /3 mode are strongly enhanced 
by matter effects. 

For sin^ 2^13 <0.1 and neutrino trajectories crossing the core of the earth the transition 
probability z/^ <->• z/e as the function of energy has three peaks: two MSW peaks due to the 
resonance enhancement of the oscillations in the mantle and in the core and the parametric 
resonance peak between them. The parametric peak dominates over the MSW peaks for 
sin^ 26*13 < 0.06. 

For sin^ 26*13 > 0-15 the energy intervals for the MSW resonances in the core and mantle 
strongly overlap and complex interference phenomena occur. 

The resonance effects manifest themselves in the zenith angle dependences of the charged 
leptons produced in the interactions of neutrinos, although the integration over the neutrino 
angles and energies leads to a smearing of the zenith angle distribution of charged leptons. 

We have found that the asymmetry A^/'^(0.6,1) of the multi-GeV e-like events can be 
as large as about 0.28 for the domain of parameters allowed by the GHOOZ bound and 
S 23 = 0.75. This value is achieved at Am ^2 — 10”^ eV^ and sin^ 2^13 > 0.6. The parametric 
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resonance is not operative at such large values of sin^ 26 * 13 . 

The asymmetry in the parameter region sin^ 26*13 ~ 0.06, Am ‘^2 ~ (1 — 2) x 10“^ eV^ 
where the parametric resonance becomes important (contributes > 60%) can be as large as 
A^/^( 0 . 6 ,l) ~ 0 . 22 , i.e. is close to the maximal possible value. 

An important consequence of the parametric effects is that even for very small values of 
the mixing angle 6*13 quite a sizeable asymmetry of the multi-GeV e-like events can result. 
The asymmetry can be as large as 0.15 even for sin^ 26*13 — 0.02. 

The oscillations of the sub-GeV neutrinos could be parametrically enhanced only for 
about a factor of 3 smaller values of Amgg, close to the lower bound of the range allowed 
by the SK data. However even in this case the parametric resonance would not affect the 
asymmetries of the sub-GeV data signihcantly. 

We have found that taking into account the subdominant <->■ oscillations leads 
to an excess of e-likes events and improves the £t of both multi-GeV and sub-GeV e-like 
data in the SK experiment. However, for all allowed values of the oscillation parameters the 
predicted excess (~ 3 — 5%) is smaller than the observed one. Thus, if the observed excess 
survives future experimental tests, one will need alternative explanations for it. 

The following remarks are in order. 

1. If the excess of multi-GeV e-like events is at least partly due to the z/g ^ oscilla¬ 
tions, it leads to a lower bound on the mixing angle 6*23: 

s 23 >i/f~ 0 . 4 , or 023 >39°. (37) 


This is an independent conhrmation of the conclusion that z/^-z/t- mixing should be rather 
large which follows from the data hts performed by the SK collaboration [lOT. However, 
unlike in the SK analysis, the lower bound (l37|) does not depend on the p-like data and is 
therefore complementary. If the excess of the sub-GeV e-like events is also at least partly 
due to the z/g z/^ oscillations, an even more stringent limit on 6*23 would follow: 6*23 ^ 43°. 

2. With increasing statistics of the SK experiment and new independent measurements 
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of the primary cosmic ray flux it will become possible to give a definitive answer to the 
question of whether there is a non-vanishing up-down asymmetry of the multi-GeV electron 
events beyond the small asymmetry due to the geomagnetic effects. If the answer is positive, 
this would be a signature of non-zero mixing 613 . In a signihcant domain of the allowed 
values of the asymmetry is strongly enhanced by the resonance matter effects on 

neutrino oscillations and in particular, by the parametric resonance. This makes it possible 
to probe even very small values of ^13 in the atmospheric neutrino experiments. 

The determination of the value of ^13 from the up-down asymmetry of the e-like events 
will require an independent knowledge of the values of Aml 2 and 6 * 23 . These could be 
obtained, e.g., from various samples of the /i-like events which only weakly depend on 6 * 13 . 

3. In principle, it is possible to experimentally disentangle the contributions from dif¬ 
ferent resonance structures. Although this does not seem to be possible with the presently 
available data, such an analysis may become possible with future data with better statistics 
and/or more accurate reconstruction of neutrino energies and directions. 

Clearly, selecting events in the non-vertical bins in which the trajectories that do not 
cross the earth’s core dominate will allow one to estimate the effects of the MSW resonance in 
the mantle. With high statistics e-like data this can be a realistic task. Clear identihcation of 
the MSW effect on oscillations of atmospheric neutrinos in the earth would be of paramount 
importance. 


For vertical bins one can use various energy cuts to discriminate between different res¬ 
onance effects. Indeed, as we have pointed out in Sec. 3, the SK response function for 
multi-CeV events has a steep low-energy slope. Therefore using different energy cuts for 
leptons one can exclude the effects of the low energy peaks in the oscillation probability. For 
instance, if Am^ ~ 2 x 10 “^ eV^, then the response function corresponding to the threshold 
1.33 CeV will cover all three resonances. However, with the threshold of 2 CeV, the effect 
of the parametric peak can be signihcantly suppressed. 


4. It is interesting to note that if the preliminary BESS results [13] are conhrmed and 
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the overall normalization of the atmospheric z/g and fluxes is indeed somewhat below the 
current theoretical predictions, the SK data would imply a smaller dehciency of atmospheric 
z/^’s. At the same time, the data on e-like events as well as the observed small value of the 
double ratio R would therefore mean a signihcant excess of the e-like events. In such a 
situation the excess of e-like events due to the parametric effect can be enhanced. 

Indeed, a smaller dehciency of atmospheric z/^’s means that the value of 6 ^ 23 , though still 
in a range of large mixing, may be farther away from 45°. If ^23 is noticeably larger than 
45°, the excess of e-like events is further enhanced by the — 1) factor [see (p!4D]. Thus, 
in this case of the reduced hux normalization, the parametric effects would be especially 
important for understanding the SK atmospheric neutrino data. 

Note added : 


After this work had been practically accomplished, the paper has appeared in which 
the SK data were analyzed in the three neutrino oscillation scheme. The zenith angle 
distributions are shown for large ^ 13 , where parametric resonance effects play no role. The 
£t of the data in ||^ agrees with our results: In fig. 16 of that paper, the allowed region in 
the triangle for Am ‘^2 = 10“^ eV^ corresponds to S 23 > 0.5 and small nonzero ^ 13 . 
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Appendix : Cross-sections and event rates 


To describe the neutrino cross section we have considered separately the processes of quasi¬ 
elastic scattering, single pion production and multi-particle production |2^. We have also 
included nuclear effects according to the treatment of Smith and Moniz [^] - the nucle¬ 
ons bound in the oxygen nucleus were assumed to fill a Fermi sphere up to a maximum 
momentum pf = 220 MeV, and to have a binding energy of 25 MeV. 


The quasi elastic scattering was described following Llewellyn Smith [^, using Fa^Q"^) = 
— 1.25 (1 -|- Q'^/Ma)~'^ for the the axial-vector form factor, with Ma = 1.0 GeV |2^. The 
nuclear effects are important for the quasi-elastic cross section: the processes where the hnal 
state nucleon is scattered in an occupied state are prohibited by the Pauli blocking effect 
and the cross section is reduced. The Fermi momentum of the bound nucleons also has the 
effect of broadening the angular distribution of the hnal state charged leptons. 


The cross section for the single pion production in the region W < 1.4 GeV {W is the 
mass of the hadronic system in the hnal state) was described following Fogli and Nardulli 


26| . In this region the most important dynamical ehect is the presence of the A resonance. 

All the other scattering processes were described using the standard formula for deep 
inelastic scattering using the leading order parton distribution functions (PDF’s) of Gluck, 
Reya and Vogt In the monte carlo calculation we used the LUND algorithms p8|, to 
construct physical particles from the hadronic state composed of the scattered (anti-)quarks 


and the nucleon remnants were described as a gg, gggg or gggg system. 


We have used a monte carlo method for the calculation. This allows to include all the 
dynamical features in detail, including the important nuclear ehects, and also to simulate 
(at least crudely) the experimental selection criteria (in particular the ‘single-ring’ and 
containment conditions). Our monte carlo also generated neutral-current events, but we 
have not considered the possibility of the mis-classification of NG events as GG events. 


To simulate single-ring events we selected the events with a charged lepton in the appro- 
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priate range of momentum, and required in addition the absence of photons or additional 
charged particles above the Cherenkov threshold. The single ring requirement is important 
because it preferentially selects lower energy neutrinos and therefore changes the response 
function for the different categories of events. For the containment requirement, we assumed 
that all electron events in the hducial volume were contained; for each muon event generated 
a neutrino interaction point in the hducial volume and checked whether the range in water 
of the hnal state /i^ was shorter or longer than the distance along the trajectory from the 
interaction point to the PMT surface. 

Our no-oscillation calculation is approximately 20% lower in normalization than the 
Super-Kamiokande monte carlo for all hve categories of events (e-like and /i-like sub-GeV 
and multi-GeV fully contained events and partially contained events), with very good agree¬ 
ment in the zenith angle distributions. The absolute normalization of the calculation is 
sensitive to details such as the minimum amount of Gherenkov light that a charged parti¬ 
cle needs to have in order to produce an additional visible ring. For the purposes of our 
discussion we hnd the agreement to be good. 
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Figure captions 


Fig. 1. Dependence of the transition probability P 2 on neutrino energy for ^m \2 = 2 x 10“^ 
eV^ (a) for sin^ 26*13 = 0.025 and different values of cos©,^ : —0.98 (solid curve), —0.88 
(long-dashed curve), and —0.85 (short-dashed curve); (b) for cos©,^ = —0.98 and different 
values of sin^ 26 * 13 : 0.014 (short-dashed curve), 0.025 (long-dashed curve), 0.057 (solid curve). 
Neutrino energy in eV. 

Fig. 2. Zenith angle dependences of multi-GeV events, (a) e-like events: solid line - no 
oscillations, dashed line corresponds to oscillations with sin^ 26*13 = 0.06, sin^ 6*23 = 0.75 
and Am |2 = 10 “^ eV^. (b) /i-like events: upper solid histogram is for no-oscillations case, 
the dashed histogram - oscillations with sin^ 26*13 = 0.06, sin^ 6*23 = 0.75 and Am '^2 — 10“^ 
eV^, the lower solid histogram - two-neutrino oscillations with the same parameters but 
sin^ 26*13 = 0. 

Fig. 3. Dependence of the up-down asymmetries of multi-GeV (a) e-like events, A^/^(0.6,1), 
and (b) /i-like events, A^/'^(0.6,1), on Aml 2 for S 23 = 0.75 and different values of sin^ 26 * 13 : 
sin^ 26*13 = 0.03 (dashed curve), 0.06 (solid curve), 0.10 (dot - dashed) curve. The dotted 
curve in £g. 3b shows the asymmetry for pure oscillations with S 23 = 0.75. (c). 

The same as in fig. 3a but for different bins: A^/^(0.84,1) (dashed curve), A^/^(0.60, 0.84) 
(dot-dashed curve) and A^/^(0.60,1) (solid curve); sin^ 2013 = 0.06. 

Fig. 4. Dependence of the up-down asymmetries of multi-GeV (a) e-like events, A^/^(0.6,1), 
and (b) /i-like events, il^/'^(0.6,1), on sin^ 20 i 3 for S 23 = 0.75 and different values of Aml 2 - 
0.5 X 10“^ eV^ (dot-dashed curve), 1.7 x 10“^ eV^ (dashed curve), 3 x 10“^ eV^ (solid curve). 
The squares on the curves represent the GHOOZ bound: parts of the curves on the right of 
the squares are excluded. 

Fig. 5. Dependence of the up-down asymmetries (a) of e-like events, 1), and (b) 

of /i-like events, A^/^(0.6,1), on S 23 for Am |2 = 1-7 x 10“^ eV^ and different values of 
sin^ 20 i 3 : sin^ 20 i 3 = 0.03 (dashed curve), 0.06 (solid curve), 0.10 (dot - dashed curve). The 
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dotted curve in fig. 5b shows the asymmetry for pure Vr oscillations with S 23 = 0.75. 

Fig. 6 . Zenith angle dependence of the multi-GeV e-like events. The solid histogram is 
for the no-oscillation case. The dashed histogram is calculated for = 1.7 x 10“^ eV^, 
sin^ 26^13 = 0.10 and sin^ 26^23 = 0.75. The points are the 535 days data of the Super- 
Kamiokande. 

Fig. 7. Iso-asymmetry contour plot for multi-GeV e-like events in the (sin^ 26 ^ 13 , Am 32 ) 
plane for S 23 = 0.75. The closed curve corresponds to = 0.264. The other curves 

(from bottom upward): = 0.15, 0.175, 0.20, 0.225, 0.25, 0.275, 0.30, 0.325, 0.35, 

0.375, 0.40, 0.425 and 0.45. The shaded area shows the region excluded by GHOOZ. 

Fig. 8 . Zenith angle distribution of the e-like events in the sub-GeV range. Solid histogram - 
without oscillations; dashed histogram - oscillations with Am ‘^2 — 0-3 x 10 “^ eV^, sin^ 2^13 = 
0.1 and sin^ 6*23 = 0.75. The points are the 535 days data of the Super-Kamiokande. 

Fig. 9. (a) The ratio of the e-like events rates in the sub-GeV region with and without 
oscillations as the function of Am ^2 ^23 = 0.75 (solid curves) and sin^ 623 = 0.65 

(dashed curves) and different values of sin^ 20 i 3 . From the lowest to the highest curve: 
sin^ 26x3 = 0.03, 0.06, 0.10, 0.20. (b) the same as in fig. 7a but for /i-like events, sin^ 26*13 = 
0.03,0.06,0.10,0.20 (from the highest to the lowest curve). The dotted curve corresponds 
to pure Ur oscillations with sin^ 6*23 = 0.5. 
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